Field dependent quasiparticles in models of strongly correlated electrons 
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In earlier work we showed how the low energy behavior of the symmetric Anderson model in a 
magnetic field H could be described in terms of field dependent renormalized quasiparticles. Here 
we extend the approach to the non-symmetric Anderson impurity model and to the infinite dimen- 
sional Hubbard model within the DMFT approach. We present NRG results for the local spectral 
densities and the local longitudinal and transverse dynamic spin susceptibilities for different param- 
eter regimes and a sequence of values of the magnetic field. We calculate renormalized parameters 
which characterize the quasiparticle excitations. Away from half-filling we find quasiparticle weights, 
z a (H), which differ according to the spin type a =f or a =j. Using the renormalized perturbation 
theory, we show that the low energy features in the NRG results can be well described in the metal- 
lic phase in terms of these field dependent parameters. We also check Luttinger's theorem for the 
Hubbard model and find it to be satisfied in all parameter regimes and for all values of the magnetic 
field. 

PACS numbers: 72.15.Qm, 75.20Hr, 73. 21. La 



I. INTRODUCTION 



A feature of strongly correlated electron systems, such 
as heavy fermions, is their sensitivity to an applied mag- 
netic field, which makes a magnetic field a useful exper- 
imental probe of strong correlation behavior. A man- 
ifestation of this sensitivity is the very large paramag- 
netic susceptibility observed in these systems. In terms 
of Fermi liquid theory, the large paramagnetic suscep- 
tibility can be interpreted as due to quasiparticles with 
exceptionally large effective masses. These large effective 
masses arise from the scattering of the electrons with the 
enhanced spin fluctuations induced by the strong local 
Coulomb interactions. An applied magnetic field sup- 
presses the spin fluctuations causing a reduction in the 
effective masses, which can be seen experimentally in de 
Haas-van Alphen measurements^. Not only do the ef- 
fective masses depend on the magnetic field, they may 
also differ for the spin up and spin down electrons^. 
Another feature that reflects the enhanced sensitivity to 
an applied field is metamagnetic behavior, where the spin 
susceptibility in a finite field H increases with the 

field strength such that dx(H)/dH > 0, which has been 
observed in some heavy fermion compounds 5 . As a re- 
lated feature it has also been predicted that strong mag- 
netic fields can induce localization in narrow conduction 
bands 6 ^, and this has been observed experimentally in 
quasi-two dimensional organic conductors 9 .. 

In an earlier paper—, we looked at local models 
of strong correlation behavior in a magnetic field and 
showed how the quasiparticles can be described by field- 
dependent renormalized parameters. We based our ap- 
proach on the single impurity Anderson model^, which 
can be characterized by the three independent param- 
eters, Ed, the impurity level, A, the broadening of this 
level due to the hybridization with conduction electrons, 
and U, the interaction at the impurity site. In the ab- 
sence of a magnetic field, the low energy quasiparticle 



excitations can be described by an effective version of 
the same model with three corresponding renormalized 
parameters, id, A, and j 7 12 i 13 i 14 . We have shown that, 
in the presence of a magnetic field H, these parameters 
can be taken as field-dependent, id, a {H), A a (H), and 
U(H). For the particle-hole symmetric model, we calcu- 
lated the TJ-dependence explicitly from numerical renor- 
malization group (NRG) calculations of the low energy 
excitations^!!. The process of de-renormalization of the 
quasiparticles with increase of mangetic field can be fol- 
lowed in these results. The parameters which describe 
the quasiparticles change slowly on increasing the field 
from zero and revert to their uncorrelated mean field val- 
ues in the extreme high field limit. 

The renormalized parameters are not just a convenient 
way of describing the low energy behavior; they com- 
pletely specify the model. A renormalized perturbation 
theory (RPT) can be set up in which the free propagators 
correspond to fully dressed quasiparticle s 15 ' 16 . This for- 
malism is particularly effective for describing the Fermi 
liquid regime, as only diagrams up to second order have 
to be taken into account to obtain asymptotically exact 
results for the T = susceptibilities, and the leading 
T 2 term in the conductivity. This perturbation expan- 
sion is not restricted to the low energy and low tem- 
perature regime, and can be used for calculations on all 
energy scales. We have shown that a very good descrip- 
tion of the T = spin and charge dynamics for the An- 
derson model in the Kondo regime can be obtained by 
summing the RPT diagrams for repeated quasiparticle 
scattering^. The results give an accurate description of 
the spin and charge susceptibilities for arbitrary mag- 
netic field values H, and for frequencies to extending over 
a range significantly larger than the Kondo temperature 
Tic- The Korringa-Shiba relation— and the sum rules for 
the spectral density are satisfied. 

In this paper we further extend this renormalized pa- 
rameter approach to infinite dimensional lattice models, 
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concentrating particularly on the Hubbard model. How- 
ever, before doing so, we generalize some of our earlier re- 
sults for the Anderson model to situations away particle- 
hole symmetry. We do this not only for the sake of com- 
pleteness but also so we can make a comparison with the 
results for the lattice model in this regime. 



II. THE NON-SYMMETRIC ANDERSON 
MODEL IN A MAGNETIC FIELD 

The Hamiltonian for the Anderson models is 

Ham = ^ £ d,o-4A + Un d ^n dil (1) 

a 



k.a 



where ed.a = £d — is the energy of the local- 

ized level at an impurity site in a magnetic field H, U 
the interaction at this local site, and T4 jCr the hybridiza- 
tion matrix element to a band of conduction electrons 
of spin a with energy e k . a — c9cPbH/2, where g c is the 
g-factor for the conduction electrons. When U = the 
local level broadens into a resonance, corresponding to 
a localized quasi-bound state, whose width depends on 
the quantity A a (uj) = TrJ2 k \V k ^\ 2 S(uj - e k ,a)- For the 
impurity model, where we are interested in universal fea- 
tures, it is usual to take a wide conduction band with a 
flat density of states so that A a (w) becomes independent 
of u>, and can be taken as a constant A a . In this wide 
band limit A a {oj) will be independent of the magnetic 
field on the conduction electrons, so we can effectively 
put g c = 0. When this is the case A CT is usually taken to 
be a constant A independent of a. 

In the renormalized perturbation theory approac h 15 i 16 
we cast the corresponding Lagrangian for this model 
£am(£c1,<t, A, U) into the form, 

^AM^d.cr, A, U) = £AM(£d,<r, A a , U) + £ c t(Al, A 2 , A 3 ), 

(2) 

where the renormalized parameters, id. a- and A a . are de- 
fined in terms of the self-energy S CT (a;) of the one-electron 
Green function for the impurity state, 



G„(w) 



1 



(3) 



lu - eaa + iA - £ CT (u;) ' 

and are given by 

ed,<7 = 2a(ed,<7 + £<r(0)), A a — z a A 1 (4) 

where z a is given by z a = 1/(1 — £^.(0)). The renormal- 
ized or quasiparticle interaction U, is defined in terms 
of the local total 4-vertex T<h(l>i,u}2, ^3>k>4) at zero fre- 
quency, 




-4 -2 
ln(h/7iA) 

FIG. 1: Plots of the renormalized parameters, Af(/i)/A, 
A|(/i)/A, s d (h)/nA, U(h)/U, fj(h), for the asymmetric An- 
derson model, with irA — 0.1, U/ttA — 2 and Ed/irA — —0.3, 
as a function of the logarithm of the magnetic field h/irA. 
The ratio A|(/i)/Af (h) is also shown. 



It will be convenient to rewrite the spin dependent quasi- 
particle energies in the form, id, a — id{h) — ahrj(h), 
where 



(6) 



where i d (h) and f)(h) are both even functions of the mag- 
netic field h — gnBH/2. 

The renormalized perturbation expansion is in powers 
of the renormalized interaction U for the complete La- 
grangian defined in equation @. The counter term part 
of the Lagrangian C ct (\i,\ 2 , A3) essentially takes care of 
any overcounting. The parameters, id,a, A, T and U, have 
been taken to be the fully renormalized ones, and the 
counter term parameters, Ai, A2 and A3, are required to 
cancel any further renormalization. They are completely 
determined by this condition. 

We can calculate the renormalized parameters by iden- 
tifying the Lagrangian £AM(£d,a, A a , U) on the right 
hand side of equation @ with the leading corrections 
to the fixed point in a numerical renormalization group 
calculation for the same modelii^ii&S ,. Results for the 
general model in the absence of a magnetic fielcU^, and 
for the symmetric model in the presence of a field, were 
given earlier—. There are some significant differences for 
the case of a magnetic field in situations without particle- 
hole symmetry which will describe here briefly. The main 
difference is that the wavefunction renormalization fac- 
tor z a (h) depends on the spin index a, and as a conse- 
quence so does the effective resonance width A a (h) 1 so 
the equations given earlier for the particle-hole symmet- 
ric model 10 have to be generalized. The induced magne- 
tization M(h) is given by M(h) = <?/iBm(/i), where 



L> = z T z i r Ti (0, 0,0,0). 



(5) 



m{h) = i(n dT -n dl ) = ^ cr tan 1 ) > !,i 
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which can be derived from the Friedel sum rule. The 
longitudinal susceptibility Xi(M ( m units of (<?/!b) 2 ) is 
given by 

Xi{h) = 0.25(p T (0, h) + p t (0, h) + U{h)pi(0, ft)p;(0, h)), 

where p a {ui,K) is the free quasiparticle density of states 
given by 



p a {u,h) = - 



(9) 



The corresponding transverse susceptibility Xt{h) (zero 
applied field limit in the transverse direction) is given by 



Xt(h) 



m{h) 
2h 



(10) 



The total occupation of the impurity site n(h) — (n^-t + 
Tidi) can be derived similarly, and is given by 



n{h) 



1 



-E 



tan 



(11) 



2tt^ \A„(h) 

and the local charge susceptibility Xc(h) is given by 

Xc(h) = 0.25[p T (0, h) + ft(0, h) - U(h)ft(0, fc)pi(0, fc)]. 

(12) 

In figure 1 we display the renormalized parameters as 
a function of the magnetic field on a log scale for the 
bare parameters Ed/^A = —0.3 and U/irA = 2, corre- 
sponding to a impurity occupation in the absence of a 
field, (n<i,a) — n(0)/2 ~ 0.35. The overall trend is very 
similar to that for the particle-hole symmetric case in 
the strong coupling regime. We do see, however, that 
Af(ft) ^ A;(/i), except asymptotically as h — > and 
h — » oo. We note that, though Aj(/i) increases mono- 
tonically with increase of h, A^(/i) initially decreases. In 
this case, where the impurity level is less than half-filled, 
the ratio Af (h)/ A i (h) > 1. This ratio is reversed, so 
A|(/i)/A|(/i) < 1, when the impurity level is more than 
half-filled. Shown in figure 2 is the result for m(h) de- 
rived by substituting these parameters into equation ([7]) , 
compared with results obtained by the direct evaluation 
of the d-site occupation values in the ground state as de- 
termined from the NRG. There is a small but systematic 
difference, of the order of 2%, between the two sets of 
results. The corresponding estimates of the average oc- 
cupation number n(h)/2 as a function of magnetic field 
h are shown in the same figure. For this quantity the 
two sets of results are almost indistinguishable. In the 
extreme large field limit the average occupation of the 
impurity level tends to unity, as the majority spin level 
gets pulled further and further below the Fermi level, 
and the impurity becomes completely polarized. In this 
regime the average renormalized level, id, approaches the 
mean field value, £d = £d + 0.5Un(h), rather than the 
bare value s&, while the other renormalized quantities 
approach their bare values. 




-4 

ln(h/7lA) 

FIG. 2: The induced magnetization m(h) as a function of the 
logarithm of the magnetic field h for the asymmetric Anderson 
model with the same set of parameters as given in figure 1. 
The dashed curve is that calculated from the direct evaluation 
of the occupation values from the NRG ground state, and the 
full curve is that deduced from the renormalized parameters 
in equation (0. Also shown is the average occupation n(h)/2 
as calculated from the NRG ground state (dashed curve) and 
the quasiparticle occupation values (full curve) as given in 
equation {TTJ. 
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FIG. 3: The longitudinal and transverse impurity site sus- 
ceptibilities, Xi(h) and Xt(h), as a function of the logarithm 
of the magnetic field h for the asymmetric Anderson model 
with the same set of parameters as given in figure 1. Xi(h) is 
calculated from equation (|8]) and Xt(h) = m(h)/2h. 



The longitudinal and transverse spin susceptibilities, 
Xi{h) and Xt(h), are plotted in figure [3] as a function of 
the logarithm of the magnetic field h. xz(M is calculated 
from equation ([5]) and Xt(h) = m(h)/2h as calculated 
from equation ([?])■ They should asymptotically converge 
to the same result in the limit h — > 0. There seems to be 
a very small discrepancy, of about 1%, between the two 
estimates in this limit. 

The generalization of our earlier results^ for the dy- 
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FIG. 4: The imaginary part of the longitudinal dynamic sus- 
ceptibility, Xi( w > ft) at T = as a function of ui, for the asym- 
metric Anderson model for h/nA = 0.001, and the same set 
of parameters as given in figure 1. The dashed curve is calcu- 
lated from a direct NRG calculation and the full curve from 
equation (|13p with the renormalized parameters. 
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FIG. 5: The spectral density for the spin up (majority) elec- 
trons pi(u>, h) (full curve) and spin down (minority) electrons 
Pl(ui,h) (dashed curve) at T — as a function of u for the 
same set parameters and magnetic field value as used for the 
plot of Xi( LJ i h) in figure 4. The two curves are almost coinci- 
dent on the scale shown. 



namic longitudinal susceptibility Xl{ u i taking account 
of repeated quasiparticle-quasihole scattering, is 

Xi{u,h)= (13) 
Xtf{u,h) + xii(v,h) + 4Ui(h)xii{u,h)xu(ui,h) 
2(l-4U?(h)x u (u,h)x u (u,h)) 

where the analytic expression for x aa (w,K) is given in 
the Appendix. The zero frequency irreducible particle- 
hole vertex Ui(h) in this scattering channel is given by 



Ui(h) = 



i + J[i + TP{h)[fr + PiW uHh){~ P] - p{f 



2(f>T+p l + U(h)p 1 p l ) 



where we have simplified the notation, /5 CT (0, h) = p a . In 
the absence of a magnetic field, or with a magnetic field 
for the particle-hole symmetric model, P] — p[, the result 
simplifies to Ui(h) = U(h)/[1 + U(h)p(0, h)], which is the 
value used in the earlier workii. 

The imaginary part of the longitudinal dynamic sus- 
ceptibility is shown in figured] for a magnetic field value 
h/irA = 0.001 for the same set of parameters, as for fig- 
ures 2 and 3. The dashed curve is that from a direct 
NRG evaluatio n 21 ' 22 and the full curve is that calculated 
using equation (|13p . with the corresponding renormal- 
ized parameters. In the direct NRG evaluation we have 
used the improved metho d 23 ' 24 of evaluating the response 
functions with the complete Anders- Schiller basis 2 ^ so 
that the sum rule for the total spectral density is satis- 
fied exactly. There is very good agreement between the 
two sets of results. The Korringa-Shiba relation does not 
hold for the model without particle-hole symmetry in the 
presence of a magnetic field. Where it does hold, in the 
absence of a field, or with particle-hole symmetry, the 



renormalized perturbation expression satisfies it exactly. 
The spectral densities, p-\{u>) and Pi(u>), as calculated 
from the NRG for this field value and the same set of 
parameters, are shown in figure [5] For this value of the 
field there are only a small differences between the two 
spectral densities; they are almost coincident on the scale 
shown. 

The dynamic transverse susceptibility Xti^i h), taking 
into account the quasiparticle-quasihole repeated scatter- 
ing, is given by 



Xt(w,h) 



xn( w > h ) 



i-u t (h)xn(u,h)' 



(15) 



where the analytic expression for x\i(w, h) is given in the 
Appendix. The zero frequency irreducible particle-hole 
vertex Ut(h) in this channel is 



U t (h) 



4/> 2 + (A T - A x ) 



m(h) 



(A T -Aj.) , ( (h-e d ) 2 +A 2 r 
2nm(h) 111 I (h+e d ) 2 +Al 



Ah 



(16) 

where h = hrj(h), and m(h) is given in equation l[7]). 
This expression simplifies in the case, Aj = A^, to give 
Ut(h) — (h — h)/m(h), which is the same as that used 
earlier—. 

Apart from a factor of 2, the dynamic transverse sus- 
ceptibility in weak fields is similar in form to the longi- 
tudinal one. In stronger fields, however, one of the peaks 
in the imaginary part is suppressed while the other peak 
is enhanced. Results are shown for the imaginary part of 
Xt(w, h) in figure[6]for a field value h/irA = 0.04 and the 
same parameters as for figures 2-5. 

The full curve gives the results derived using equation 
(fT5)l . and the dashed curve from a direct NRG calcu- 
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FIG. 6: The imaginary part of the transverse dynamic sus- 
ceptibility, h) at T = as a function of w, for the asym- 
metric Anderson model for h/nA = 0.04 and the same set of 
parameters as given in figure 1. 



FIG. 8: The spectral density for the spin up electrons (uj, h) 
(full curve) and spin down p[(u),h) (dashed curve) at T — 
as a function of u> for the same set parameters and magnetic 
field value as used for the plot of x*( w > h) in figure 7. 
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FIG. 7: The imaginary part of the longitudinal dynamic sus- 
ceptibility, Xi( w > ft) at T = as a function of ui, for the asym- 
metric Anderson model for the same parameters and magnetic 
field value (h/nA = 0.04) as given in figure 6. 

lation. The peak positions are in good agreement, and 
the slightly broader peak from the NRG data can be at- 
tributed to the logarithmic broadening used in the di- 
rect NRG evaluation. There is a sum rule, that the total 
spectral weight is equal to —2m(h), which is satisfied pre- 
cisely both in the RPT result and also in the NRG cal- 
culation, as we have used the improved prescription for 
the response functions based on the complete Anders- 
Schiller basis. The imaginary part of the correspond- 
ing dynamic longitudinal susceptibility \i(uj,h) for this 
value of the magnetic field is shown in figure [7J Due 
to the stronger magnetic field, the peaks are somewhat 
suppressed as compared with the results shown in fig- 
ure for h/irA = 0.001, but the overall features are very 
similar. The spectral densities, pi(u>) and Pi(w), for this 
field (h/nA = 0.04) are shown in figure [51 As the po- 
larization is much stronger for this higher magnetic field 



value, there is now a marked difference between the two 
spectral densities for the two spin types. 

In summary, the results of this section for the impurity 
models away from particle-hole symmetry display new 
features, such as the spin dependent resonance widths 
in the quasiparticle density of states. We show that the 
magnetisation and the static response functions can be 
well described in terms of the field dependent renormal- 
ized parameters. We have used these parameters to cal- 
culate the dynamic longitudinal and transverse spin sus- 
ceptibilities in the RPT formulae and we find excellent 
results when compared with those obtained from a direct 
NRG calculation. 



III. INFINITE DIMENSIONAL HUBBARD 
MODEL 

In this section we turn our attention to lattice mod- 
els with a local interaction terms and in particular to 
the infinite dimensional Hubbard model. The Hubbard 
model has played a similar role for lattice models as the 
Anderson model for impurity models, being the simplest 
model of its type, where the interplay of kinetic energy 
and strong local interactions can be studied. 



A. Dynamical mean field approach and 
renormalized parameters 

The Hamiltonian for the model in a magnetic field is 
given by 

H = - /XUjclyCj^ +h.c.) - y] ahn ia + 11}^ n»,t n »,.l» 

icr i 

(17) 
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where Uj are the hopping matrix elements between sites 
i and j, U is the on-site interaction. 

From Dyson's equation, the Green's function Gk )ff (w) 
can be expressed in the form, 



Gk,«r(w) 



1 



-£ CT (k,^)-e(k)' 



(18) 



where So-(k, uj) is the proper self-energy, and e(k) = 
e~ k '( Ri ~ Rj )tjj , \i a = ji + cr/i, with the chemical po- 
tential /i. The simplification that occurs for the model 
in the infinite dimensional limit is that So-(k, uj) becomes 
a function of to onl y 26 ' 27 . In this case the local Green's 
function G 1 ° c (uj) can be expressed in the form, 



G l ° c (u) = ^G k ,,rH = J de 



D(e) 



U + H a - S CT (cj) - £ 



(19) 

where D(e) is the density of states for the non- interacting 
model (U = 0). It is possible convert this lattice problem 
into an effective impurity oneSS and write this Green's 
function in the form, 



where 



(20) 



(21) 



The Green's function G 1 ° c (uj) can be identified with the 
Green's function G a {uj) of an effective Anderson model, 
by re-expressing Qq^uj) as 



so that 



Q 1(uj) = uj + /i + ah - K a ((j), 



G a {u>) = 



OJ — £da - K a (uj) - £ ff (u>) : 



(22) 



(23) 



with Sda = — Mo-- The function K a (uj) plays the role of 
a dynamical mean field describing the effective medium 
surrounding the impurity. In the impurity case in the 
wide band limit we have K a (uS) = — iA. Here as can be 
seen from equations ([2T1l and ([20|) . K a (uj) is a function 
of the self-energy T, a (u>), and hence depends on a. As 
this self-energy is identified with the impurity self-energy, 
which will depend on the form taken for K a (uj), it is clear 
that this quantity has to be calculated self-consistcntly. 
Starting from an initial form for K a (ui), Y, a (ui) is calcu- 
lated using an appropriate 'impurity solver' from which 
G 1 ° c (uj) can be calculated using equation (fT9|) . and a new 
result for K a (uj) from equations (f2Tj) and (|22|) . This 
K a {oS) serves as an input for the effective impurity prob- 
lem and the process is repeated until it converges to give 
a self-consistent solution. These equations constitute the 
dynamic mean field theory (DMFT), and further details 
can be found in the review article of Georges et ali^. 



Of the many impurity solvers the most commonly used 
are the Monte Carlo, the Exact diagonalization (ED) 
method and the NRG, all of which have advantages and 
disadvantages. Here, we wish to calculate the field de- 
pendent renormalized parameters to describe the quasi- 
particles, as we did for the Anderson model, so we use the 
NRG approach. It is also the most accurate method for 
calculations at T = and for the low energy excitations. 
There has been a DMFT study of the static properties of 
a half-filled Hubbard model in a magnetic field using the 
ED method by Laloux et alA The focus of our paper 
here, however, is rather different so there is little over- 
lap with this earlier work but, where there is, we make 
comment and compare with their results. 

We need to specify the density of states D(uS) of the 
non-interacting infinite dimensional model, which is usu- 
ally taken to be either for a tight-binding hypercubic or 
Bethe lattice. Here we take the semi-elliptical form cor- 
responding to a Bethe lattice, 



ttD 



(24) 



where 2D is the band width, with D = 2t for the Hub- 
bard model, and /^o the chemical potential of the free 
electrons. We choose this form, rather than the Gaus- 
sian density of states of the hypercubic lattice, as it has 
a finite bandwidth. 

In the NRG approach 2 ^ the conduction band is log- 
arithmically discretized and the model then converted 
into the form of a one dimensional tight binding chain, 
coupled via an effective hybridization V a to the impurity 
at one end. In this representation K a {uj) — |V r CT | 2 go,<r( w )i 
where <7o,<r(w) is the one-electron Green's function for the 
first site of the isolated conduction electron chain. If we 
use this form in equation (|23|) . and expand the self-energy 
E CT (o;) to first order in uj, we can write the result in the 
form, 



G a (uj) = 



UJ — Sdo 



|^| 2 5 o, CT M + 0( W 2 ) 



2V 



(25) 



where 



e Aa = Zcr {e Aa + £<r(0)), \V a \ 2 = ZvlVvl 2 , (26) 

where z a = 1/(1 — £^.(0)), very similar to the result in 
equation ^ . We can interpret this result as a free quasi- 
particle propagator, Gq >( j{uj), given by 



Go,«t(w) 



uj - e A a- \V a \ 2 g ^{uj) 



(27) 



and the quasiparticle weight. A similar Fermi liquid 
expansion in (fT9|) leads to 



(28) 



J uj + ^ 0: <T - e 



where /2o.c = z (7 (fi a — S CT (0)). In the DMFT approach 
we identify 



(29) 
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which specifies the form of g ^ a (u) in (|27|) and and yields 
/2o,ct = — £da- We can define a density of states po tCr (u>) f° r 
the free quasiparticles via po ja (u>) = — ImG?o, -(a> + i<5)/7r. 
For the Bethe lattice, this takes the form of a band with 
renormalized parameters, 



irDl 



L>1 - (w + /x , CT ) 2 



(30) 



We can then compare the results based on these RPT 
formulae, which take into account the repeated quasi- 
particle scattering, with the NRG results over the whole 
frequency range. 

Having covered the basic theory, we are now in a po- 
sition to survey the results for the Hubbard model in 
different parameter regimes. 



where D a = z a D. 

We can deduce the renormalized parameters ida and 
V a from the low lying excitations calculation from the 
the NRG using a generalization of the method given in 
an earlier paper—. The quasiparticle weight z a is then 
obtained from the relation z a = |V^/V^-| in equation 

(HI. 

We can define a quasiparticle occupation number n° by 
integrating this density of states up to the Fermi level, 



ft 



dw po,<rM- 



(31) 



Using Luttinger's theorem 2 - it is possible to relate this 
free quasiparticle occupation number n° to the expecta- 
tion value of the occupation number n a in the interacting 
system at T = 0. Using the quasiparticle density of states 
in equation (|30p . we can rewrite equation (|3ip as 



dwD(w)8(/i ff -E 5 (0)-w), 



(32) 



where 6(u>) is the Heaviside step function and D(uS) as 
given in equation J24|) . Then from Luttinger's result the 
right-hand side of equation (f32|) is equal to n a . We then 
have the result, 



(33) 



that the occupation for electrons of spin a is equal to the 
number of free quasiparticle of spin a, as calculated from 
equation (|3"Tj) . It should be noted that there is no simple 
generalization of the h = result, p, — jiQ = 2(0), in the 
spin polarized case. 

We can also calculate the local longitudinal and trans- 
verse dynamic spin susceptibilities, Xi( w ) an d Xt{u), for 
this model using equations, (|13p and (|T5|) . derived earlier 
in the impurity case. The free quasiparticle susceptibili- 
ties Xa,cr' ( w ) required are calculated using the local free 
quasiparticle density of states po ja (u) given in equation 
([30)1 . We can calculate the local on-site quasiparticle in- 
teraction U as in the impurity case, but we do not have 
the simple formula relating U to x;(0) and Xt(0) that 
enabled us to deduce the irreducible quasiparticle inter- 
actions Ui and Ut ; the impurity formula we used earlier is 
only valid in the wide band limit. To determine Ui and Ut 
in the lattice case we use the condition that Re%; (u>) and 
Rext(w) fit the NRG result at the single point u> = 0. 



B. Results at Half-filling 

We present the results at half-filling for three main 
parameters regimes where we find qualitatively different 
behavior. The results in all cases will be for a Bethe 
lattice with a band width W = 2D = 4, setting t = 1. 
In concentrating on the field induced polarization, we do 
not include the possibility of antiferromagnetic ordering. 
The regimes are a relatively weak coupling regime (a) 
where U is smaller than the band width, an intermediate 
coupling regime (b) with W < U < U c , where U c is 
the value at which a Mott-Hubbard gap develops in the 
absence of a magnetic field (U c ~ 5.88), and (c) a strong 
coupling regime with U > U c . 

The first plot in figure [5] gives the spectral densities for 
the majority spin electrons p-\(u>) for various magnetic 
field values in the weakly correlated regime, U — 2. We 
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FIG. 9: The local spectral density for the majority spin Pi(ui) 
for U — 2 and various fields h. 

can see clearly that, for increasing magnetic field, more 
and more spectral weight is shifted to lower energies (the 
opposite happens for the other spin component, which 
is not displayed here). Above h ~ 1.0 the system is 
completely polarized, 2m = 1. This extreme high field 
limit corresponds to an insulator; there is a gap of the 
magnitude A s (/i) = 2h + U — IF between the upper (mi- 
nority) and lower (majority) bands, which both have the 
semi-elliptical form as for non-interacting system with 
W = 4. The inverse of the quasiparticle weight z a (h), 
which corresponds to the enhancement of the effective 
mass m*(h) = m/z ff (h), is shown as a function of h in 
fig. [TUl The renormalized parameters (RP) £d,<r an d \V<j\ 2 
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were calculated from the NRG low energy excitations, as 
described in the previous section, and the values of z a (h) 
deduced using equation (|26|) . We also calculated z a (h) 
from the numerical derivative of the NRG calculated self- 
energy S ct (oj) for both spin up and spin down electrons. 
All values agree very well and z-\(h) = Z[(h), as expected 
in the case with particle-hole symmetry. The method 
based on the renormalized parameters is only applica- 
ble in the metallic regime. The values of z a (h) increase 
from about 0.75 to 1, which corresponds to a progressive 
"de-renormalization" of the quasiparticles with increas- 
ing field, as observed earlier for the impurity models. 
Since the interaction term in the Hubbard model acts 
only for opposite spins it is clear that there is no renor- 
malization when the system is completely polarized with 
one band fully occupied and the other empty. The ex- 
pectation value of the double occupancy (n^n{) decreases 
with increasing field, which further illustrates why the in- 
teraction term becomes less important for larger fields. 




° 5 0.2 OA 0.6 0.8 1 

h 

FIG. 11: The renormalized chemical potential fio, a (h) calcu- 
lated from renormalized parameters (RP) and directly from 
the self-energy for U = 2 and various fields h. 
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FIG. 10: The inverse of the quasiparticle weight z a (h) calcu- 
lated from renormalized parameters (RP) and directly from 
the self-energy and the magnetization m(h) for U = 2 and 
various fields h. 

We can also follow the field dependence of the renor- 
malized chemical potential //q )0 -(/i) as shown in figure QT] 
ft is shown deduced from the renormalized parameter 
(RP) £d,o- and as calculated directly from the self-energy. 
The agreement is very good over the full range of mag- 
netic fields. 

Mean field theory is valid for very weak interactions, 
so we compare our results for /2o,<t(M for U = 2, with 
the mean field value = \i + ah — Un™ l a in figure [11] 
The results coincide for h = 0, when /i™^ = and when 
the system becomes fully polarized at large field values, 
/lo 1 ^ = —a(U + h), but in general ji™ f a > fio,<j(h). We also 
compare the mean field (MF) result for the field depen- 
dence of the magnetization m(h) with the one obtained 
in the DMFT calculation in figure [T2j The general be- 
havior is similar, but the mean field theory without quan- 
tum fluctuations overestimates the magnetization, as one 
would expect. 




X) 0.2 0.4 0.6 0.8 1 



m(h) 

FIG. 12: The magnetization in the mean field approximation 
compared with the DMFT result for U = 2 and for the full 
range of magnetic fields h. 



In the next plot in figure [T3l where U = 5, we show 
typical behavior in the intermediate coupling regime. 
Similar to the weak coupling regime, we find a shift of 
spectral weight towards lower energy for the majority 
spin. There is, however, a difference in the way this hap- 
pens due to the initial three peak structure, namely the 
quasiparticle peak in the middle gets narrower for in- 
creasing field and finally vanishes in the polarized phase. 
The quasiparticle weight, which is shown in figure [T4l re- 
flects this behavior by decreasing to zero with increasing 
field signaling heavy quasiparticles. When the material 
is polarized, however, the z a (h) reverts to 1, which corre- 
sponds to the band insulator as before. This approach to 
the fully polarized localized state in high fields contrasts 
with that found in the weak coupling regime. It gives 
rise to metamagnetic behavior in this parameter regime. 
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FIG. 13: The local spectral density for the majority spin p-\{u)) 
for (7 = 5 and various fields h. 
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FIG. 15: The real part of the local longitudinal dynamic spin 
susceptibility for (7 = 5 and various fields h. 
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FIG. 14: The inverse of the quasiparticle weight z a (h) calcu- 
lated from renormalized parameters (RP) and directly from 
the self-energy for (7 = 5 and various fields h. The inset shows 
the magnetization m(h). 



To illustrate further this different response to a mag- 
netic field, the real part of the local longitudinal dynamic 
spin susceptibility Xli^i h) as a function of w is shown for 
various values of h. It can be seen that the local suscepti- 
bility x loc (/i) = Kexi(0, h) in this regime increases with 
h so that dx loc {h)/dh > 0. This can also be seen in 
the curvature of the magnetization shown in the inset of 
figure [TJl This is behavior characteristic of a metamag- 
netic transition and related to the magnetic field induced 
metal-insulator transition. 

We can also check the Luttinger theorem in a magnetic 
field, as discussed in the previous section, by comparing 
the values of n£ , deduced from integrating the quasiparti- 
cle density of states with the value of n a calculated from 
the direct NRG evaluation in the ground state. These 
results are shown in figure [TBI It can be seen that there 
is excellent agreement between the results of these two 
different calculations, ffi a — n CT , so that Luttinger's theo- 
rem is satisfied for all values of the magnetic field in this 



intermediate coupling regime. 




0.05 0.1 0.15 0.2 



h 

FIG. 16: The comparison of the spin dependent occupation 
numbers n° and n„ corresponding to Luttinger's theorem in 
a magnetic field for (7 = 5 and the range of fields h. 

Finally we consider the strong coupling regime with 
U > U c , where for h = the spectral density has a Mott- 
Hubbard gap so that for half-filling the system is an insu- 
lator. The electrons will be localized with free magnetic 
moments coupled by an effective antiferromagnetic ex- 
change J ~ t 2 /U. In fields such that h > J, the system 
polarizes completely as can be seen in figure 1X71 where we 
show the total density of states p{ui) = P]{u) +p±(ui) for 
h = and h = 0.2. For smaller fields, such that h < J, 
we do not find a convergent solution to the DMFT equa- 
tions, and the iterations oscillate between local states 
which are either completely full or empty. We interpret 
this as due to the tendency to antiferromagnetic order 
which in a weak field, due to the absence of anisotropy, 
will be almost perpendicular to the applied field in the x- 
y plane with a slight canting of the spins in the z-direction 
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FIG. 17: The total local spectral density p(u>) for {7 = 6 for 
h — (dashed line), Mott insulator, and h — 0.2 (full line), 
fully polarised band insulator. 
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FIG. 19: The local spectral density for the minority spin pi (w) 
for U = 5, x = 0.5 and various fields h. 



(spin flopped phase) . In this calculation no allowance has 
been made for this type of ordering, but this state can 
be well described using an effective Heisenberg model for 
the localized moments. 



C. Quarter Filled Case 

We now compare the results in the intermediate cou- 
pling regime with U — 5 at half-filling with those at 
quarter filling, x = 0.5, where the Fermi level falls in the 
lower Hubbard peak in the spectral density. To see how 
the band changes with increasing magnetic field we plot 
the density of states for both spin types, for the major- 
ity spin electrons in figure [18] and for the minority spin 
electrons in figure HH for various values of the magnetic 
field. 







— h =0 


i fi 




— h =0.1 


* tl 




--h =0.25 


§ t If 
' ' // 




— h=0.4 


> // 

if // 
i'll 
1 * // 


\ 


ii / / 

■ ' / / 

■ / / 
ii / / 


il 

in 
H 

Ml 

'Ii 




1 1 / / 


Ml 




i* / / 


i u 




/ 


■ \\ 

i 1 \\ 




hi/ 
*// 


i \\ 














. . . -—7 .i .... , .«t 



-2 i 4 " " 6 8 



03 

FIG. 18: The local spectral density for the majority spin Pj(w) 
for U = 5, x — 0.5 and various fields h. 

In the majority spin case the lower peak broadens on 
the low energy side and the weight in the upper peaks 



decreases with increase of the field. The opposite features 
can be seen in the minority spin case, with the spectral 
weight in the lower peak below the Fermi level decreasing 
and the weight in the upper peak increasing. Thus the 
increase of spectral weight below the Fermi level for the 
majority spin electrons, and the decrease for the minority 
spin electrons, can be seen to be due to a change of band 
shape rather than a simple relative shift of the two bands, 
which would be the case in mean field theory. In the 
fully polarized state there are no minority states below 
the Fermi level and the upper peak in the majority state 
density of states has disappeared. 

The corresponding values for the inverse of the quasi- 
particle weight 1/ z a (h) are shown in figure [^01 for a range 
of fields. As noted in the impurity case, the quasiparticle 
weights differ for the two spin types with Zj(/i) > zi(h). 
The values of z a (h) have been calculated, as described 
earlier, both from the energy levels (RP) and from a nu- 
merical derivative NRG derived self-energy. There is rea- 
sonable agreement between the two sets of results, and 
the small differences to be seen be attributed to the un- 
certainty in the numerical derivative of the NRG self- 
energy. As in the impurity case, there is an initial de- 
crease of zi(h) with increase of h, whereas z^(h) increases 
monotonically. The field dependence of the magnetiza- 
tion is also shown in figure [20l and is similar to the half- 
filled case with a weak interaction (U = 2). We have 
calculated, but do not show, the corresponding occupa- 
tion values for n° which again agree well with the values 
of iv, verifying Luttinger's theorem. 

Our conclusion from these results, and from calcula- 
tions with other values of U , is that when there is sig- 
nificant doping, the behavior in the field corresponds to 
a weakly correlated Fermi liquid, very similar to that at 
half-filling in the weak interaction regime. The only re- 
markable difference in the field is the spin dependence 
of the effective masses as shown in figure 1201 which was 
already found similarly in the impurity case. 
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FIG. 20: The inverse of the quasiparticle weight z a (h) calcu- 
lated from renormalized parameters (RP) and directly from 
the self-energy for U = 5, x — 0.5 and various fields h. The 
magnetization m(h) is also displayed. 



D. Near half filling 

Very close to half-filling and for large values of U we 
have a qualitatively different parameter regime. Here the 
system is metallic but we can expect strong correlation 
effects when U is of the order or greater than U c , due 
to the much reduced phase space for quasiparticle scat- 
tering. We look at the case with 5% hole doping from 
half-filling and a value U = 6, which is just greater than 
the critical value for the metal-insulator transition. We 
show the spectral density of states for both the majority 
and minority spins states and various values of the mag- 
netic field in figures [21] and [22j respectively. There is 
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FIG. 21: The local spectral density for the majority spin p-\(u)) 
for U = 6, x = 0.95 and various fields h. 

a clear sharp quasiparticle peak for h = at the Fermi 
level at the top of the lower Hubbard band. As in the 
quarter filling case with U — 5 we see a similar trans- 
fer of spectral weight with increasing field to below the 
Fermi level for the majority spin case, and above the 
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FIG. 22: The local spectral density for the minority spin pi (w) 
for U = 6, x — 0.95 and various fields h. 



Fermi level for the minority spins. For large fields when 
the system is completely polarized there is still a sharp 
narrow peak in the spectral density of the minority spin 
states above the Fermi level, though the spectrum for 
the majority states below the Fermi level is that of the 
non-interacting system. A spin up electron added above 
the Fermi level feels no interaction as the system is com- 
pletely spin up polarized so these electrons see the non- 
interacting density of states. On the other a spin down 
electron above the Fermi level interacts strongly with the 
sea of up spin electrons. The self-energy due to scatter- 
ing with particle-hole pairs in the sea creates a distinct 
resonance in the down spin density of states just above 
the Fermi level. Just such a resonance was predicted by 
Hertz and Edwards^ for a Hubbard model in a strong 
ferromagnetic (fully polarized) state. 

The field dependence of the inverse of the quasiparticle 
weight is presented in figure 1231 Again we find reason- 




0.25 

FIG. 23: The inverse of the quasiparticle weight z a (h) calcu- 
lated from renormalized parameters (RP) and directly from 
the self-energy for U = 6, x = 0.95 and various fields h. The 
inset shows the magnetization m(h). 
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able agreement between the two methods of calculation 
for these quantities. The magnetization as a function of 
h is shown as an inset in the same figure. The behavior 
of Zf(h) and z^(h) as a function of h contrasts sharply 
with the behavior found for the metallic state at half- 
filling with [7 = 5 shown in figure HU For zero field 
the quasi-particle weight has a very similar value in both 
cases. At half-filling the tendency of the magnetic field to 
induce localization resulted in values of and 
{zf(h) = zi(h)) which decrease sharply as a function of 
h. In the 5% doped case with U — 6, the system must re- 
main metallic and the quasiparticles weights, z^{h) and 
Zi(h), both increase in large fields though their values 
differ significantly. The quasiparticle weight for the mi- 
nority spin electrons decreases initially with increase of h, 
whereas that for the majority spins z\ (h) increases mono- 
tonically and quite dramatically with h. When the sys- 
tem becomes fully polarized the up spin electrons become 
essentially non-interacting, z^(h) — 1, whereas there is a 
strong interaction for a down spin electron and we find in 
this case Zi(h) ~ 0.15. The interpretation for this is as 
given in the previous paragraph for the spectral densities. 

We conclude that already a small doping of the system 
is enough to maintain a metallic character even for very 
strong interaction. Although the zero field spectra of 
the half filled case for [7 = 5 and the small doping case 
with [7 = 6 display very similar zero field behavior, i.e. a 
strongly renormalized quasiparticle band with similar z a , 
no field induced localization transition occurs for finite 
doping and no metamagnetic behavior is observed in the 
latter case. 



IV. QUASIPARTICLE DYNAMICS 

Having deduced the renormalized parameters of the 
quasiparticles from the NRG results presented in the pre- 
vious section, we are now in a position to test how well 
we can describe the low energy dynamics of this model 
in a magnetic field in terms of a renormalized perturba- 
tion theory. We look at the various parameter regimes in 
turn. 



A. Free quasiparticle spectral density 

It is of interest first of all to see how the free quasiparti- 
cle density of states f>a(uS) multiplied by z a (h) compares 
with the spectral density p a (u>). In figure [2"4"1 we make a 
comparison in the zero magnetic field case. We see that 
the quasiparticle band gives a good representation of the 
low energy peak in p a {ui) and, as expected, does not re- 
produce the high energy features. These, however, to a 
fair approximation can be described by the mean field so- 
lution p m {(w) weighted with a factor 1 — z a as can be seen 
in figure |2"4"1 A case with a finite magnetic field h = 0.15, 
where the peaks in the density of states of the two spin 
species are shifted due to the induced polarization rela- 
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FIG. 24: The free quasiparticle density of states in comparison 
with interacting local spectral density for (7 = 5 and h = 0. 
We have also plotted a thin black line for p m f (w) = [D(ui + 
U/2) + D(lu — U /2)]/2 which describes the non-magnetic mean 
field solution and weighted with 1 — z„ . 



tive to the Fermi level, is shown in figure [221 The figure 
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FIG. 25: The free quasiparticle density of states in comparison 
with interacting the local spectral density for (7 = 5 and h = 
0.15. 

focuses on the region at the Fermi level and one can see 
the the free quasiparticle density of states describes well 
the form of p a (ui) in the immediate vicinity of the Fermi 
level. It is to be expected that the frequency range for 
this agreement can be extended if self-energy corrections 
are included in the quasiparticle density of states using 
the renormalized perturbation theory as shown in the 
impurity case^i. 

In the fully polarized case with h — 0.22 there is 
complete agreement between the quasiparticle density of 
states and p a (u>) for both spin types as can be seen in 
the results shown in 1261 where z-\ = z\_ = 1. 

In figure |2"T1 we show a fully polarized result [h — 0.26) 
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FIG. 26: The free quasiparticle density of states in comparison 
with interacting local spectral density for U = 5 and h = 0.22. 



for the case near half filling, x = 0.95, U = 6 discussed 
in section III.D. We can see that the different values for 




FIG. 27: The free quasiparticle density of states in comparison 
with interacting local spectral density for U — 6, x — 0.95 and 
h = 0.26. 



the field dependent quasiparticle weight for up and down 
spin Za-(h) lead to remarkably different quasiparticle band 
shapes. With Zf ~ 1 the majority spin quasiparticle band 
is essentially that of the non-interacting density of states. 
The very much smaller value z± leads to a narrow quasi- 
particle band above the Fermi level. The low energy flank 
of this quasiparticle band describes well the narrow peak 
seen in the spectral density just above the Fermi level. 
To describe these strong asymmetries in the spectral den- 
sities near half filling, we need z-f 3> Zj,, which contrasts 
with the cases at half filling such as in figures [25] and [26] 
where always Zj — z \- This suggests a discontinuous be- 
havior of the renormalisation factors z a on the approach 
to half filling. 



B. Dynamic susceptibilities at Half-filling 

We now compare the NRG results for the longitudi- 
nal and transverse local dynamic spin susceptibilities for 
the same value U = 5 and a similar range of magnetic 
field values with those based on the RPT formulae (fl"3"]) 
and (fl5|) . In figure [28] we show the imaginary part of the 
transverse spin susceptibility calculated with the two dif- 
ferent methods. It can be seen that RPT formula gives 




FIG. 28: A comparison of the imaginary parts of the trans- 
verse dynamic spin susceptibility for U = 5 and h — 0.0 
calculated using the renormalized parameters (RP) and from 
a direct NRG calculation. 

the overall form of the NRG results, and precisely fits 
the gradient of the NRG curve at u> — 0. Some of the 
relatively small differences between the results might be 
attributed to the broadening factor used in the NRG re- 
sults which gives a slower fall off with to in the higher 
frequency range, and a slightly reduced peak. We get 
similar good agreement between the two sets of results 
for the same quantity for the case with a magnetic field 
h = 0.15, shown in figure [29] 

In figure 1301 where we give both the real and imagi- 
nary parts of the transverse susceptibility for h — 0.19, 
we see that this overall agreement is maintained in the 
large field regime where we get metamagnetic behavior. 
The shapes of the low energy peaks for both quantities 
are well reproduced by the RPT formulae. Note that the 
peak in the real part is not at w = 0, so it is not fixed 
by the condition that determines Ut, but nevertheless is 
in good agreement with the NRG results. Due to their 
very small values it becomes difficult to calculate z a (K) 
as the system approaches localization for larger fields. In 
this regime as z a {h) — > the free quasiparticle density 
of states will converge to a delta-function. Self-energy 
corrections to the free quasiparticle propagators, which 
were used in the calculation of will become in- 

creasing important as this limit is approached. Once the 
system has localized, and is completely polarized, how- 
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FIG. 29: Plots of the imaginary part of the transverse dy- 
namic spin susceptibility for U = 5 and h — 0.15. 



FIG. 31: The real part of the longitudinal dynamic spin sus- 
ceptibility for U = 5 and h — and h — 0.15. 
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FIG. 30: The real and imaginary parts of the transverse dy- 
namic spin susceptibility for U — 5 and h — 0.19. 
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FIG. 32: The imaginary part of the longitudinal dynamic spin 
susceptibility for (7 = 5 and h = 0.15. 



ever, we saw in figure [2H that the values deduced of Ji a 
(z a (h) = 1) gave the a quasiparticle density of states 
coinciding with the NRG result. 

Results for the longitudinal susceptibility are shown in 
figures [31] and [32j In figure [31] we give the values for 
the real part as a function of ui for h = and h = 0.15. 
Here the peak height, which is at uj — 0, is fixed by the 
condition which determines Ui . The width of the peaks in 
the two sets of NRG results, however, is given reasonably 
well by the RPT equations. The imaginary part of the 
longitudinal susceptibility obtained by the two methods 
is given in figure l3"2l for h = 0.15. Again there is overall 
agreement between the two sets of results. The slight 
undulations seen in the RPT results are due to the sharp 
cut off in the band edges in the free quasiparticle density 
of states. For larger values of h the agreement with the 



NRG results is not as good as that as for the transverse 
susceptibility, and the central peak in the real part of the 
RPT results narrows more rapidly with h than in those 
obtained from the direct NRG calculation. 



C. Dynamic susceptibilities at quarter filling 

We give two examples of results for the susceptibilities 
for the model at quarter filling for the case U = 5 and h = 
0.1. In figure [33] we give the real and imaginary parts for 
the transverse susceptibility. Despite the large value of U, 
we can see that the peak heights are very much reduced 
compared with those seen in the half- filled case for U = 5. 
The peak widths are also an order of magnitude larger 
as can be seen from the w-scale. Nevertheless the RPT 
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FIG. 33: The real and imaginary parts of the transverse dy- 
namic spin susceptibility for (7 = 5, x = 0.5 and h = 0.1. 



results reproduce well the overall features to be seen in 
the NRG results. The real and imaginary parts for the 
longitudinal susceptibility are shown in figure [331 
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FIG. 34: The real and imaginary parts of the longitudinal 
dynamic spin susceptibility for U = 5, a; = 0.5 and h — 0.1. 

Again all the low energy features are reproduced in the 
RPT results. In this regime, apart from the overall factor 
of 2, there is less difference between the transverse and 
longitudinal susceptibilities than at half-filling. 
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FIG. 35: The real and imaginary parts of the transverse dy- 
namic spin susceptibility for U = 6, x = 0.95 and h = 0.15. 



seen on an w-scale an order of magnitude smaller than 
that for quarter filling due to the much stronger renor- 
malization effects in this regime. There is excellent agree- 
ment both with the peak positions and shapes between 
the NRG and RPT results for both quantities. This is 
also seen to be the case for the real and imaginary part of 
the longitudinal susceptibility shown in figure [551 though 
the peak in the real part can be seen to be marginally 
narrower in the RPT results. 
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FIG. 36: The real and imaginary parts of the longitudinal 
dynamic spin susceptibility for U = 6, x = 0.95 and h = 0.15. 



D. Dynamic susceptibilities near half filling 

Near half-filling (x = 0.95) we show plots for the two 
susceptibilities for parameters U = 6 and h = 0.15. In 
figure [35] we give the real and imaginary parts of the 



V. SUMMARY 

In this paper we have extended our earlier work-i^, 
where we described the low energy behavior of the sym- 
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metric Anderson model in a magnetic field h in terms of 
field-dependent renormalized quasiparticles, to the non- 
symmetric Anderson model. The main new feature that 
emerges is the dependence of the quasiparticle peak res- 
onance width A a (h) on the spin type a as well as on the 
value of the magnetic field. The T = spin and charge 
susceptibilities can be expressed as explicit formulae in 
terms of the parameters A a (h), £<i, a (h) and a local field 
dependent interaction between the quasiparticles U(h). 
It was also shown earlier for the symmetric modelii, that 
an excellent description of the low energy spin dynamics 
can be obtained by taking into account repeated quasi- 
particle scattering in a renormalized perturbation theory 
(RPT). We have shown that these results can also be gen- 
eralized to the non-symmetric model for the transverse 
and longitudinal spin susceptibilities, which again agree 
remarkably well with those obtained from a direct NRG 
calculation. 

The method for calculating the field dependent quasi- 
particle parameters has also been extended here to in- 
finite dimensional lattice models where the self-energy, 
as in the impurity case, is a function of frequency only. 
We have applied these methods to the Hubbard model 
to describe the low energy excitations in terms of a spin 
dependent renormalized quasiparticle band. 

For the Hubbard model at half filling, where z-\{h) = 
z±(h), we presented results for the three main parameters 
regimes where the model displays qualitatively different 
behavior. Our results are on the whole consistent with 
those obtained earlier by Laloux et al. 6 obtained using 
the ED method. An exception is in the insulating regime 
for weak fields, where we could not find a convergent 
solution of the DMFT equations. We attributed this to 
the fact that in this regime the magnetic field is smaller 
than the exchange coupling between the localized spins 
so that ground state would be one in which the spins 
would have a canted antiferromagnetic ordering in the 
plane perpendicular to the field. 

Well away from half filling we find a magnetic response 
similar to the weakly correlated case even for large values 
of U. The large phase space for quasiparticle scattering in 
this regime leads to modest renormalization effects. Here, 
as in the impurity case, we find spin dependent quasipar- 
ticle weights, &f(h) ^ z±(h). This implies spin depen- 
dent as well as field dependent effective masses, which 
have been discussed earlier in work by Spalek et 
and Riseborough^ 3 .. The calculations by Spalek et al. 
were based on a Gutzwiller 32 and a mean-field slave bo- 
son approach^. We can make a comparison of our results 
(section III.D) near half filling, x = 0.95, with theirs in 
the later worki. We find a qualitatively similar behavior 
with the majority spin effective mass decreasing with h, 
but quantitatively there are differences. The field depen- 
dence of the minority spin effective mass l/z±(h) shows 
a very slow increase initially in both sets of results, but 
the large field behavior is quite different. As seen in fig- 
ure [53] we find a significant decrease in l/zj(/i) for large 
fields whereas the corresponding quantity in figure 3 in 



reference^ increases. 

The strong magnetic field dependence of the effective 
masses found in the calculations by Riseborough is based 
on the assumption that the system is close to a ferromag- 
netic transition (paramagnon theory). However, DMFT 
calculations for the Hubbard model find that any ferro- 
magnetism in the Hubbard model only occurs in a very 
small region of the parameter space near half-filling and 
for very large values of UM- Our results are well away 
from this regime and the large effective masses obtained 
here can be attributed to the tendency to localization 
rather than the tendency to ferromagnetism. 

Using the field dependent renormalized parameters 
z a (h) and jX a {h) in the RPT formulae for the dynamic 
local longitudinal and transverse spin susceptibilities we 
found agreement with the overall features to be seen in 
the NRG results for these quantities. In the case of 
the transverse spin susceptibility excellent agreement was 
found in all the metallic regimes and for all values of the 
magnetic field considered, except in the high field regime 
at half filling as the localization point is approached, 
where consistent values of the renormalized parameters 
are difficult to calculate. The comparison of the RPT 
results with those from NRG was also excellent for the 
longitudinal dynamic susceptibility in the weaker field 
regime h < 0.15 but less good for higher fields, h > 0.15. 

In all metallic parameter regimes a spin dependent 
Luttingcr theorem in the form n a = n°, the number of 
particles equals the number of quasiparticles, was found 
to be satisfied for all strengths of the magnetic field. In 
this form it even holds in the fully polarized insulating 
state. 

Phenomena like field and spin dependent effec- 
tive masses and metamagnetic behavior have been 
observed experimentally in several heavy fermion 
compounds^ ' 3 ' 5 ' 35 . The Hubbard model, however, being 
a one band model is not an appropriate starting point to 
make a quantitative comparison with the heavy fermion 
class of materials. A periodic Anderson model with a two 
band structure and including the degeneracy of the f elec- 
trons would be a better model to describe these materials. 
Field dependent effects in this model have been studied 
by several techniques, modified perturbation theory 3 ^, 
exact diagonalization^l, 1/JV expansion^ and variational 
approach 39 . The approach used here could be generalized 
to the periodic Anderson model, but restricted to the 
non-degenerate case and N = 2 as it is computationally 
too demanding in the NRG to deal with higher degen- 
eracy. The Hubbard model at half filling has been used 
as a lattice model to describe the strongly renormalized 
Fermi liquid 3 He£i However, the metamagnetic behav- 
ior predicted for relevant parameter regime is not seen 
experimentally^. In section III.D we found for small 
doping large effective masses, but no metamagnetic be- 
havior. This raises the possibility that the weakly doped 
Hubbard model could serve as a basis for interpreting the 
experimental results for liquid 3 He. 
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VI. APPENDIX 

The free quasiparticle dynamic susceptibility X<r,<r' ( w ) 
for the impurity model in the wide band limit, Aj = 
Aj, were given earlier—. Here we give the more general 
results for Aj ^ Aj, 



i 



ttuj u - 2iA„ 



.=-1,1 \ 



aid,a + «A CT 



(34) 



for lo > 0, and for u = 0, 

W(0) = /5 a (0). 



(35) 



The values for u> < follow from the fact 
that Rex^M = ReXo-^-w) and Imx^.o-M = 
-Imx<7,o-(— w). For a' ^ a, 
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